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DISTRIBUTIONS  WITH  APPLICATIONS 
by 

S.  Iyengar  and  Y.L.  Tong 


1.  Introduction. 

The  problem  of  evaluating  multivariate  probabilities  arises  in  many 
areas  of  statistics,  such  as  in  the  construction  of  confidence  regions  for 
the  mean  vector  or  regression  parameters  in  a  general  linear  model,  the 
determination  of  critical  regions  for  certain  tests,  multiple  comparisons, 
reliability,  etc.  Because  these  probabilities  are  hard  to  evaluate  numeri¬ 
cally,  a  thorough  study  of  the  nature  of,  say,  the  distribution  function  is 
needed  to  suggest  good  approximations  and  inequalities.  For  example, 
tables  for  selected  values  of  parameters  are  available,  and  one  usually 
interpolates  linearly  to  get  an  approximation;  if  we  know  that  the  tabulated 
function  is  convex  or  concave,  we  also  get  a  bound  for  the  quantity  of 
interest . 

It  is  well  known  that  the  cumulative  and  rectangular  probabilities  for 
elliptically  contoured  distributions  are  increasing  functions  of  the  correla¬ 
tions  (see  Tong  [8]  for  more  precise  statements),  but  the  rate  of  increase 
has  not  yet  been  studied.  In  this  paper,  we  prove  some  results  concerning 
the  convexity  in  correlations  of  the  distribution  functions  of  random  vectors 
with  elliptically  contoured  distributions,  and  their  absolute  values.  Sections 
2  and  3  concentrate  on  the  bivariate  case,  and  Sections  4  and  5  provide  some 
applications.  The  last  section  contains  extensions  to  the  multivariate  case; 
the  results  here  are  not  as  strong  as  for  the  bivariate  case,  as  the  computa¬ 
tions  are  much  harder. 

•  1 


2.  The  Bivariate  Normal. 

The  following  identity  due  to  Plackett  [6]  facilitates  the  computa¬ 
tions  that  follow. 


Lemma  1.  Let 
and  variances 

(1) 


(f>n(x;$)  be  the  multivariate  normal  density  with  mean  zero 

i.  if  I =  (pij)»  then 

3d)  =  32<t> 

3pij  '  ' 


When  n»*2,  we  write  <£(x;P)  instead  of  <J>2(x;$).  Identity  (1)  gives 
rise  to  natural  conditions  on  partial  correlations  for  certain  results  to 
be  true.  Now  let 


(2)  F2(x;p)  *  P(XX  <  xrX2  <  x2) 

and 

G2(x;p)  =  P(|xi|<x1,|x2|<x2) 


where  X  =  (X^,X2)'  has  density  4>(x;p).  Then  we  have 

Theorem  1.  F2  is  a  convex  (concave)  function  of  p  in  an  interval  I 

2 

and  only  if  h(x;p)  =  (px^-x^)  (px^x^)  +  p  (1-p  )_>(<)  0  holds  for  pel 
G2  is  a  convex  (concave)  function  of  pel  if  and  only  if 

2pXlX2 

2~ 

k(x;p)  =  (px1-hc2)  (px2~x1)  -  p  (1-p2)  +  e  lp  (p  (l-p2)+(px1~x2)  (px^x^ ) 


for  pel. 


•  (  <  )  0  , 


if 


2 


Proof.  Using  (1)  and  doing  the  following  interchange  of  differentiation 


and  integration,  we  get 


-x-  v 
3  p  2 


37^  *(t!0)dt 


<f>(x;p)  , 


So  that 


9p2  c2 


- 5-5-  4»(x;p)h(x;p)  . 

(1-P  ) 


A  similar  computation  yields 


3  2 

— r  G,  =  - t-t-  <})(x  ,-x_;p)k(x;p)  . 

ap  u-pV 

Now  the  theorem  is  immediate. 

More  helpful,  perhaps,  is  the  following: 

Corollary  1.  Let 


and 


If  x^x^  >  0,  F2  is  convex  in  at  least  (0,m);  if  x^x2  <  F2 
concave  in  at  least  (M,0).  G2  is  convex  in  at  least  (-m,m). 

In  particular,  if  c  =*  x^  =  x2  then  F2  is  convex  in  (0,1)  and  when 

c  >  0,  G2  is  convex  in  (-1,1).  For  F2  we  can  get  even  more  detailed 
information:  if  |c|  ^>/2-l,  F2  is  convex  in  (-1,1)  since  in  that  case 
h  is  positive  for  all  p.  In  the  next  section,  we  will  generalize  the  above 


3 


theorem  and  corollary  to  arbitrary  elliptically  contoured  distributions, 
but  the  detailed  information  in  this  paragraph  is  much  harder  to  get. 

Corollary  1  shows  that  the  interval  of  convexity  depends  on  how 
spread  out  x^  and  x ^  are.  A  more  precise  statement  of  this  uses 
the  concept  of  majorization;  for  a  full  account,  see  Marshall  and  Olkin 
[5]. 

2 

Theorem  2.  Let  x,y  eH  .  If  F^Cxjp)  is  a  convex  (concave)  function  of 
pel,  and  if  x  )>  y  (y  ^  x) ,  then  ^(yjp)  is  also  convex  (concave)  in 
I. 


Proof .  Suppose  that  F^(x;p)  >_  0  and  x  ^  y.  Note  that  h(x;p)  in 
Theorem  1  is  a  Schur  concave  function  of  x.  Thus  h(y;p)  h(x;p)  0, 

g2  ^ 

and  — y  F?(y;p)  =  - —r  (y ; p)h(y ; p)  0.  The  other  case  is  similar. 

3p^  1  (1-pV 

Corollary  1  also  shows  that  the  convexity  of  depends  upon  the 

spread  of  x^  and  X2*  Indeed,  we  have  ample  numerical  evidence  to  show 
that  Theorem  2  also  applies  to  G2;  however,  we  do  not  yet  have  an  analytical 
proof  of  it,  as  the  algebra  is  quite  cumbersome. 


3.  Bivariate  Elliptically  Contoured  Distributions. 

Some  of  the  above  results  can  be  extended  to  all  elliptically  contoured 
distributions.  Recall  that  a  random  vector  X  eRn  is  elliptically  contoured 
if  there  is  a  y  for  which  (X-y)  has  characteristic  function  f(t)  =  <j)(tf|t) 
where  (jr.tO,0^)  -+  R;  when  X  has  a  density,  it  is  of  the  form 
III  1x),  where  x11  ^g(x^)dx  <  °°.  A  representation  theorem  of 

Cambanis,  et  al.  [1]  says  that  this  is  equivalent  to  saying  that 
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(3) 


X  =  y  +  R^2!^ 


where  R  0  is  independent  of  is  uniformly  distributed  on 

l/2 

the  unit  sphere,  and  $  is  the  symmetric  square  root  of 

To  start  our  investigation,  we  need  a  generalization  of  Plackett's 
identity.  Joag-dev,  et  al.  [4]  proposed  the  following:  if  g  is  a 
function,  then 


(4) 


3f 


n 


d  ,r  ik 

’  V£ 


where  $  ^  *  (p^)  and  f(x;$)  *  |{|  ^^g(xf|  ^x) .  We  propose  another 
generalization:  suppose  X  has  density 


n 


f  (x;|,B)  -  2^l  1/28P+2  (X,r1x)exp(-Bx,r1x)  , 
H »  P  r«  /H  i  _  \ 


(5) 


n,p 


snr(f +p) 


where  s  is  the  surface  area  of  the  unit  ball  in  Kn.  Then,  we  have  the 
n 

following  lemma: 


Lemma  2.  If  X  has  density  f  ,  then  its  characteristic  function  is 

-  n,P 


<t>  (tf$t;B)  where 
n,p 


<J>  (u;  8)  =  e 

n,p 


'U/4B  !  (?)  kii1  (f  + u'V  ^)k 

k=0  i=0  Z 


and  the  density  satisfies 


(6) 


3pij  n>P  2B  m=0  m-  i=m 


js  l  Sr1"1  <f+1> 


-1 


9xi9xj  n’m 


Proof .  See  the  appendix. 
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Now  define  and  as  in  (2)  of  Section  2.  The  next  theorem 

generalizes  Corollary  1.  The  proof  consists  of  two  steps:  first,  use  (6) 
to  establish  the  theorem  for  such  densities,  and  second  use  a  geome¬ 
trical  argument  and  (3)  to  show  that  the  theorem  is  true  for  all  ellip- 
tically  contoured  distributions. 

Theorem  3.  Corollary  1  is  true*  for  all  bivariate  elliptically  contoured 
distributions. 

Proof.  Step  1.  Using  (5),  we  get 


-  P  SB  - 

dp  2  26 


m=u  *  — w  w 


m=0 


and 


Similarly,  writing  x  =  (x^,-X2)  we  have 


and 


d2  2(px1~x2)  (px^xp 


f,  (x;$,6)  + 

y  P 


2(px1+x2)  (px2+x  ) 
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It  is  well  known  that  -p-  F0  and  p  -p-  G0  are  non-negative  (Tong  [8], 

dp  Z  dp  Z 

p.  64,  ff.).  Thus,  the  theorem  follows  for  this  case  immediately  upon 


and 


2  2 


Go. 


inspection  of  the  other  terms  in  the  expression  for  — j 

dp  1  dp- 

respectively.  (The  use  of  (4)  instead  of  (5)  yields  the  same  expressions 

2  2 
for  -4f,  and 


dp 


2  2 


dp 


2  2' 


X1  x2 

Step  2.  If  we  let  m  =  —  A  —  and  assume  x-x  >  0,  then  F«  is  convex 

X£  ^1  z  z 

in  pe(0,m)  whenever  the  density  of  X  is  f-  .  Recall  that  X  =  R$^^U, 

^  >  P 
2 

where  U  is  uniform  on  the  unit  circle,  and  R  has  density 


,  /  \  8P  P  -3r 

h  „(r)  =  — r  rF  e 
P»P  P! 


1/2 

If  e^  and  e^  are  the  standard  unit  vectors,  and  b^  =  $  e^. 


then 


-OO 

,(x;p)  -  J0  P(b'U  <  b'U  <  -Z)hp>g(r)dr 


p(x,r;p)h 

Jo  p’ 


(r)dr  . 


Since  F2  is  convex  in  pe(0,m),  we  have  for  0  <_  a, 6  and  a+6  =  1, 


0  > 


(p (x,r;ap1+6p2)-ap(x,r;p1)-ep(x;r;p2) }hp  g(r)dr 


2 

for  p^,p2c(0,m)  and  all  p,B  >  0.  Now  let  p,8  -►  00  with  p/B  =  r^; 
it  is  easy  to  see  that  y(x,rg,p)  is  convex  in  pe(0,m)  for  each  r^. 
But  since  U  is  uniformly  distributed,  y  is  just  proportional  to  the 
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lengths  of  certain  arcs  on  the  circle,  so  we  now  have  the  purely  geome¬ 
trical  fact  that  this  length  is  a  convex  function  of  p.  Thus,  the 
distribution  of  R  is  irrelevant  and  we  have  that  is  convex  in 

pe(0,m)  when  x^x2  >  ®  for  aH  eliptically  contoured  distributions. 

The  proof  for  the  other  cases  in  the  same. 

The  key  to  Step  2  above  is  that  for  the  family  f  (x;i,B),  the 

n,p 

radius  vector  ||£  ^^x||  can  be  made  to  degenerate  to  any  point  Yq  0 
by  suitable  choice  of  p  and  Q.  In  particular,  this  family  of  densities 
is  not  multivariate  unimodal. 

4.  Applications. 

Two  well-known  probability  inequalities  are 

SlepianT s  inequality:  Let  X  and  Y  be  elliptically  contoured  with 
EX  =  EY  ~  0,  with  parameters  $  and  T,  respectively,  where  $  =  (p^), 
T  *  (t^)  are  correlation  matrices  and  P^_.  >_  T  .  Then  for  all 
a  «  (ax, . . .  ,an) ,  P(X±  <_  a±;i«l, . . .  ,n)  >_  P(Y±  ^  a^i-l, . . .  ,n) . 

Dunn’s  inequality:  Let  X  ^  N  (0,t).  Then  P(|x,|  <  a.;i=l . n)  > 

- -  n  I  —  i  — 

P(|xtl  1  ai)  . 

Of  particular  interest  for  the  construction  of  simultaneous  confidence 
bounds  is  the  case  a^  =  a^  ■•••■  =  a.  These  inequalities  have  been 

used  extensively  to  provide  bounds  for  the  probabilities  of  interest.  We 
shall  show  below,  however,  that  the  bounds  are  often  far  from  the  true 
probabilities;  that  is,  the  bounds  are  bad  approximations.  We  also  show  how 
convexity  results  from  Sections  2  and  3  can  yield  very  good  approximations 
(but  not  necessarily  bounds)  for  high  dimensional  probabilities. 
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Suppose  then  that  X  is  elliptically  contoured, 
and  we  want  to  approximate  P(X^  _<  a;  i=l,...,n)  and  P(  |x^  |  _<  a;i=l, .  .  .  ,n)  . 
Let  =  E^  l{X^  _<  a}  and  =  E^  so  the  above  probabilities 

are  P(Sn  =  n)  and  P(Tn»n),  respectively.  Also,  let  f  be  the  equicorre- 
lation  matrix  whose  off-diagonal  elements  are  all  equal  to  the  average  of 
the  off-diagonal  terms  of  $.  Then  we  have 

Theorem  4.  If  p..  >  0  for  all  i,1,  then  var^S  >  var  S  .  For  all 
-  ij  —  f  n-  n 

t ,  var_T  >  var  T  . 

Z  n-  J  n 

Proof.  The  variance  of  T  is 
-  n 


var i  T  »  constant  +  2  £  P(|X  |  <  a>|x. |  <  a)  *  k({p  }). 

t  n  i<  j  i  J  ij 

Now  P(|X_J  <_  a,|Xj|  a)  is  a  convex  function  of  p^  by  Theorem  3,  so 

the  Hessian  of  k  is  a  positive  definite  (diagonal)  matrix,  and  k  is  a 

convex  function.  Being  a  symmetric  function,  k  is  also  Schur-convex;  thus 

k({p  })  >  k({p}),  where  p  =  E  p  . /(!!).  The  proof  for  S  is  similar, 
ij  —  ij  l  n 


Notice  that  ElS  =  E  S  ,  so  that  S  has  the  same  location  under  t 

T  n  «  7  n  ' 


5  L  5 

+  n  J  n 


and  but  is  more  dispersed  under  $  when  p^.  >_  0.  It  is  tempting  to 

conjecture,  then,  that  certain  probability  inequalities  obtain  for  PCS^^). 


For  an  illustration,  suppose  that  a  =  0,  so  that  ExS 


is  =  E  S  »  n/2;  in 
|  n  j  n 


this  case,  we  might  conjecture  that  P±(S  =n)  >  P  (S  -n) ,  for  p  ,  .  >  0. 

\  n  —  n  ij  “ 

Indeed,  this  is  true  for  some  special  cases,  but  not  in  general.  Even  when 

the  inequality  does  not  hold,  we  can  expect  that  approximating  Pi,(Sn=n)  by 

P  (S  “n)  will  be  better  than  using  the  Slepian  bounds,  Pj.  (S  =n)  or 
f"  n  *min  n 
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are  equicorrelation  matrices  with 


Px  (S  =n) ,  where  t  and  t 
t  n  Tmax  T 

Tmav 


1  max 

parameters  max  p 


min 

^  and  min  ,  respectively.  Table  1  below  compares 


the  two  methods;  for  the  correlation  matrices  used  for  comparison,  see 
the  Appendix. 


True 

\  approx 

a) 

.216 

.214 

b) 

.134 

.135 

c) 

.243 

.241 

d) 

.129 

.126 

TABLE  1 

Slepian  upper  bound 
.272 
.228 
.183 
.008 


Slepian  lower  bound 
.164 
.106 
.273 
.316 


Of  course,  similar  comments  apply  for  a  ^  0  and  Pi.  (T^=n) .  Also, 
we  could  use  the  concavity  results  of  the  earlier  sections  to  give  simple 
sufficient  conditions  for  the  opposite  inequalities;  we  omit  the  details. 

5.  Computations  for  the  Equicorrelation  Case. 

The  above  approximations  are  most  useful  if  they  are  easily  evaluated; 
that  is,  the  approximation  be  expressible  as  a  low  dimensional  integral, 
which  then  requires  numerical  quadrature. 

For  the  normal  distribution,  the  equicorrelated  case  is  especially 
easy.  When  p  _>  0,  it  is  known  that  a  one-dimensional  integral  suffices: 

(6)  P(X,  <  a,...,X  <  a)  = 

1  —  n  — 


*(- 


-  g 


-)  <f>  (t)dt 


10 


where  <{>  and  $  are  the  standard  normal  density  and  distribution,  respec¬ 
tively.  And  for  p  <  0,  Steck  [7]  provides  a  formula  relating  it  to  the 
p  >  0  case. 

For  other  elliptically  contoured  distributions,  the  equicorrelated 
case  does  not  yield  a  one-dimensional  integral.  Even  when  p  «  0,  the 
probability  integrals  do  not  immediately  reduce  in  dimension  because  p  =  0 
implies  independence  (within  the  elliptically  contoured  family)  only  for 
the  normal  distribution.  For  one  subclass,  though,  a  two-dimensional 
representation  is  possible.  Cambanis,  et  al.  [1]  use  a  result  of  Schoenberg 
and  von  Neumann  to  show  that  if  X  has  an  elliptically  contoured  distribu¬ 
tion,  which  can  be  embedded  in  an  infinite  exchangeable  sequence,  then  X 
has  density  |£|  ^2f(x?$  ^x)  where  =  p  f°r  i  t  j  and 


(7) 


f  (u) 


[0,«) 


(2Trr2) 


-n/2 


exp(-u/2r  )dF(r) 


and  F  is  the  distribution  of  ||^  '  a||  .  Using  this  result,  we  can  apply 

Fubini’s  theorem  to  immediately  get  that 


P(X±  1  a;i=1> • • • >n) 


[0,°°) 


-1/2 


2Trr2$|  exp(- ~xf  (r2^ )  1x)dx}dF(r)  . 


The  inner  integral  is  a  Gaussian  integral,  which  can  be  evaluated  easily 
by  (6);  thus,  we  have  that  for  scale  mixtures  of  normals,  orthant  probabili¬ 
ties  and  rectangular  probabilities  are  easily  computed.  Important  special 
cases  are  the  multivariate  t  and  Cauchy  distributions  which  are  elliptically 
contoured. 
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For  other  elliptically  contoured  distributions,  the  problem  is  harder; 
it  essentially  reduces  to  that  of  computing  the  surface  area  of  the  inter¬ 
section  of  circular  patches  on  the  sphere,  since 


P(x±  1  ai;Vi) 


P(b!U  )dF(r)  . 

J[0,»)  1  r  1 


6,  Extensions  to  Multivariate  Case. 

There  are  several  ways  in  which  to  extend  our  convexity  results  to 
higher  dimensions.  For  instance,  one  might  consider  F(£;a)  =  (X^<^ a^;V^) 

and  ask  when  F  is  a  convex  function  of  A  simpler  but  instructive 

case  is  when  $  *  $  ,  an  equicorrelation  matrix,  and  a^  =  a.  Even  in 
this  case,  the  computations  are  complicated,  so  our  results  are  somewhat 
limited. 


Theorem  5.  If  X  ^  N  (0,t  ),  then  P(X,  <  a;V,)  is  a  convex  function  of 
-  n  Tp  i—  i 

p  >_  0  when  a  <_  0.  When  a  _<  1  -  /l ,  it  is  a  convex  function  of  p  in 

its  entire  domain  (-  — ^-,1). 

n-1 

Proof.  Let  (dx)  =  dx^,...,dx  and  interpret  a  0-fold  integral  as  0. 
-  n  1  n 

Now 


F(p;a)  =  P(Xi  £  a;V^) 


•a 

j  — oo 


<{).n(x;p)  (dx)n 


and 
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dF(p;a)  _  ,n. 
dp  ~  V 


<j>  (x,a,a;p) (dx) 
n  n—z 


,nv 


(2)<J>2(a,a;p) 


at(p) 


4>  (x;p')(dx)  ? 

n- z  n—  z 


where  t(p)  =  {(l-p)/(l+p)  (l+2p)}^^  and  p'  =  p/(l+2p).  Next 


>i  fat<0> 


4>_  ,(x;p’)(dx)  , 

n— z  n—Z 


+  »2(a,a;p) 


at(p) 


4>  (x;p')(dx)  } 

n— z  n— z 


By  Theorem  1,  the  first  term  is  positive  if  |a|  >  /2-1  (for  all  p)  or 
p  ^  0  (for  all  a).  Thus,  we  only  need  to  study 


_d_ 

dp 


at(p) 


4>  (x;p')(dx)  . 

n- z  n— z 


Now  recall  that  if  G(x)  ®  g(x,y)dy,  then 


G’(x) 


b  (x) 


gx(x,y)dy  +  b' (x)g(x,b(x)) 


Repeated  use  of  it  shows  that 


jd_ 

dp 


at(p) 


4>  (x;p')(dx) 

n— z  n—  z 


1  <n:2> 


rat(p) 


(l+2p) 


2  v  2 


<f>  ,(x,at(p),at(p);p')(dx)  , 

n— 4  n— 4 


+  (n-2)at'(p) 


at(p) 


<fi  , (x,at (p ) ;p ’ )  (dx)  . 
n— J  n— J 


13 


The  first  term  here  is  always  positive;  the  second  term  is  positive  only 

when  a  <_  0,  since  t1  (p)  =  -t(p){ — £  0.  Thus,  we  have  the 

1-P  1-P 

simple  sufficient  conditions  for  the  positivity  of  d^F(p  ;a)/dp^. 

Corollary  2.  If  X  is  a  scale  mixture  of  equicorrelated  normals,  then 
P(X  £  a;V^)  is  convex  in  p  _>  0  for  a  <_  0. 

Proof.  From  (7),  we  get 


P(X±  <  a;V±) 


P(Z  <  ar;V  )dF(r) 

*  [0,«)  1  1 


where  Z  ■  (Z^,...,Zn)  ^  Nn(0,^).  Since  a  <  0  =>  ar  _<  0,  the  integrand 
is  convex  for  p  _>  0  for  each  r.  Hence,  the  integral  is  also  convex. 

Note  that  this  corollary  indicates  that  the  rather  specialized  result 
for  the  normal  (a  <  1  -  /2  =>  F(p;a)  convex  for  all  p)  cannot  hold  for 
all  elliptically  contoured  distributions,  since  we  may  have  ar  >  1  -  /2 
even  though  a  <  1  -  /2 . 

We  have  only  succeeded  in  extending  Corollary  2  to  arbitrary  elliptically 
contoured  distributions  when  n=3,  by  using  the  family  of  distributions  (5). 
If  we  let 


F(p;a) 


fn  (dx)  , 

J_ oo  n»P  n 


then 


dF (p ;  a) 

dp 


(n)  — 
2B 


m=0 


n,m,p 


f  (x,a,a;$,B) (dx) 
n ,  m  n- Z 

3  7 
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where  a  ■  ^7-  11^  ^  (—+i)  See  the  "working  rules"  for  this  family 

n,m,p  m!  i=m  2 

of  distributions  in  the  Appendix;  the  notation  below  is  explained  there. 
Next,  after  some  more  algebra. 


71 '«>**>  ■  <5>  A 

dp 


(SA  ,A  f  -A  .  T  a  ,  f  ,  ) 
,  n-1  n  n,p  n,n-l  n,k,p  n,k 


(dx) 


(x,a,a;^ ,B)  n  2 

+  (S> <n22)  3  j.(p"W)an,k,p  |  fn>k<*,a,a,a,a;M><dx)n-4 

4  S  k»0  ^  * 


-  (?)  (A)  1  « 


3  48  k-o  n,k,p  n  n,k 


A  f 


(x,a,a,a;^ ,S) 


(dx) 


n-3  ' 


The  second  term  is  obviously  positive.  The  first  term  is  positive  when 
p  >  0  because 


A 

n,n,- 


1 


2p 


(l-p)(l+(n-l)p)  - 


<  0 


and 


A  .A 
n-1  n 


L  n“2  9 

- =-= - 5-  (a(l+(n-3)p)  -  p  \  x.T  >  0 

(l-p)Z(l+(n-l)p)Z  1 


The  sign  of  the  third  term  is  largely  determined  by  A^;  here 


n  (1-p) (l+(n-l)p) 


n-3 

(a (l+(n-4)p )  -  P  l  x±) 


Thus,  when  n=3,  A^  is  negative  for  a  _<  0,  and  so  Corollary  2  extends 
to  all  elliptically  contoured  distributions  (by  the  same  geometrical  argu¬ 


ment  of  Theorem  3).  For  n  _>  4,  A  <0  and  it  is  no  longer  easy  to  see 

d2  n 

where  — —  F(p;a)  is  positive. 

dpZ 
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Appendix. 

2 

(i)  Proof  of  Lemma  2.  Let  ft  ( | s |  )  =  E  exp(is'U)  be  the  characteristic 
-  n 

function  of  the  uniform  distribution  on  the  unit  sphere  in  Rn.  It  is 
known  from  the  Gaussian  case  that 


-u/2 


00  _  n-1  -r  /  2 

nn^r  U)  ~n/?  f - 7T  dr 

o  2  /2_1  r£) 


r 

J  n 


2  r"-le-r  /2u  dr 

o  n  2 n/2_1  r(-r)  un/2 


Differentiating  both  sides  with  respect  to  u  and  rearranging  some  terms, 
we  get 


f°°  0  n+2k-l  -r  /2u  /0  /0,.  ~  rv  ^  . 

^  z  2X  r  e  ,  -u/2  n/2+k  r  ,k*  _  ,  l0,w  Ni 

ft  (r  )  — - — ■  jz  -  dr  =  e  u  \  (,)  II  (n+2j  )  (-u)  . 

J0  n  T(^)2  /2_1  •  -  1 


k  .  k-i 

n 

i=0  *  j=0 


Now  recall  that  if  X  has  density  f  (x;i,B),  then  X  =  Ri^2!], 

n,  p 

2 

where  R  has  density 

,  11 

P  +  T  ,  n  - 

i  (  ^  0  2  P  +  2  1  -6r 

h  ft(r;n)  =  - —  r  e 

p’e  r(P+f) 


so  that  the  characteristic  function  of  X  is  given  by  f(t)  =  (t'itjfO, 

n,p 


where 


<t>  n(u;e) 

n,p 


ft  (r2u)h  0(r2;n)2rdr 
n  p,  p 


r  „  ,  2S,  ,r  \  2rdr 

ft  (r  )hr,  - n - 

Jg  n  P»P  u  u 


is  given  in  the  lemma. 

Next,  to  get  the  generalization  of  Plackett's  identity,  recall  that  the 
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Fourier  inversion  formula  gives 


f 

n,p 


(x;$,6) 


1 

(2TT)n 


-ix't 

e 


<f>  (t’$t;6)dt  . 

n,p 


Differentiating  under  the  integral,  we  get 


3f 


9p 


JLlE 


ij 


28 


H 


n 


-ix’t 
e  e 


t'tt 

46 


,  n 

P  +  2 


y  (p)  — 

kVk  k  n 


,  f*t  * 

(_  “46  }  dt 


n  (y+v) 
v-0 


and 


a2f 


J!i£ 


3x^3xj 


,n 


-ixf  t 
e  e 


t'U 

46 


p  k-1  ..  ,j.  k 

I  (J)  n  (f  +  v)‘l(-  ijii)  dt 
1  3  k-0  K  v=0  P 


Thus,  we  get  the  partial  differential  equation  because  the  coefficients  are 
equal: 


+  ! 


n  (f+i) 

i=0 


<k> 


m=k 


2i  (m) 

ml  'k^ 


P-1  _ 

n  (f+j) 


-1 


k-1 

n 

i=l 


,n  ,  i \  1 

(y  +  i) 


which  can  be  established  by  induction.  An  empty  product  is  to  be  inter¬ 
preted  as  1. 


ii)  For  convenience,  we  list  some  "working  rules"  for  the  family  of  densities 

f  (x;i,6).  Write 
n,p 


A 


C 

P 


2BP  2/snr(!+p)l*|1/2 


and 
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A  =  — - — a  A  =  _ _ _  A 

i  ax .  ’  ij  ax .  ax.  ’ 

i  i  j 


n 


where  s  is  the  surface  area  of  the  unit  sphere  in  H  . 


n 


a)  f  (x;$,B)  =  C  A*  e 
n,p  p 


p  -BA 


b)  f  /f  i  =  BA/(f+p-l) 
n,p  n,p-l  2 


c)  -5 —  f  =  BA.  ( - ^ -  f  f  ) 

3x±  n,p  in  +  p_1  n,p-l  n,p 


d)  -5 — 5 —  f  =  B{(BA.A  -A..)f - E -  (2BA.A.-A, , ) f  . 

9xi9xj  n,p  i  j  ij  n,p  |.+ p  _  i  1  J  U  n,p-l 


+  BA  A  - 


1  j  (|+p-l)(|+p-2)  n’P  2 


} 


e)  If,  in  Lemma  2,  we  write 


».,p  _i_ 

■  26 


„:0  n,n,,p  5xi9xj  n>" 


then 


satisfies 


p!  Pn^  ,n  .  .-1 
a  =  J-j-  II  (t+  i) 

n,m,p  m!  ±=m  2 


a  a  .  =  a  .  and  a0  =  1  . 

n,m,p  n,k,m  n,k,p  2,m,p 


iii)  The  correlation  matrices  used  to  get  the  entries  in  Table  1  are  the 
following . 
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a) 

'l.OOO  0.600  0.636  0.472' 

b) 

'1.000  0.585  0.177  0.197^ 

1.000  0.383  0.458 

1.000  0.209  0.217 

1.000  0.705 

1.000  0.291 

1.000 

1.000, 

c) 

'1.000  0.579 

0.241' 

d) 

'l.OOO  0.717 

-0.232" 

1.000 

0.582 

1.000 

-0.471 

1.000; 

1.000 

/ 
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